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1.  Introduction 

Energetic  (~10^-10^  keV)  ions  are  often  observed  in  the  magnetosheath.  Two 
dominant  sources  are  the  magnetosphere  [e.g.,  Hones  et  ai,  1972;  West  and  Buck, 
1976;  Sibeck  et  ai,  1987;  Fuselier  et  ai,  1991]  and  the  solar  wind  energized  at  the  bow 
shock  [e.g.,  West  and  Buck,  1976;  Asbridge  et  ai,  1978;  Crooker  et  ai,  1981;  Fuselier 
et  ah,  1991].  To  understand  the  origins  of  energetic  magnetosheath  ions  would  greatly 
facilitate  the  determination  of  the  mass,  momentum,  and  energy  transfer  processes 
within  the  geospace. 

The  escape  of  energetic  magnetospheric  particles  into  the  magnetosheath  at  times 
occurs  on  interconnected  magnetosphere-magnetosheath  field  lines  that  result  from 
magnetic  reconnection  at  the  magnetopause  [e.g.,  Speiser  et  ai,  1981;  Scholer  et  al, 
1981].  However,  it  mostly  occurs  on  a  continuous  basis  through  a  leakage  process  of 
finite  gyroradius  effects  [Sibeck  et  al,  1987].  This  mechanism  predicts  that  energetic  ions 
leak  from  the  magnetosphere  at  postnoon  local  times  and  energetic  electrons  at  prenoon 
local  times  due  to  the  difference  in  their  drift  paths  in  the  magnetosphere.  These  ions 
stream  along  magnetosheath  magnetic  field  lines.  Depending  on  the  magnetic  geometry 
they  may  escape  into  the  upstream  region  of  the  bow  shock  [e.g.,  Sarris  et  al,  1976; 
Luhmann  et  al,  1984]. 

The  upstream  and  downstream  regions  of  the  c|uasi-parallel  bow  shock  (the  angle 
between  the  average  upstream  magnetic  field  and  the  shock  normal,  $Bn,  less  than 
45°)  are  often  populated  with  energetic  particles  and  magnetic  waves  and  turbulence 
[e.g.,  Paschmann  et  al,  1979;  Greenstadt  at  el,  1980;  Bonifazi  and  Moreno,  1981; 
Crooker  et  al,  1981;  Ipavich  et  al,  1981;  Greenstadt,  1985;  Luhmann  et  al,  1986; 
Mobius  et  al,  1987;  Gosling  et  al,  1989;  Fuselier  et  al,  1995,  and  references  therein]. 
Beginning  with  the  earliest  observations,  it  was  realized  that  the  connection  of  the 
interplanetary  magnetic  field  (IMF)  to  the  bow  shock  is  a  necessary  condition  for  the 
presence  of  energetic  particles  on  the  field  line  [e.g.,  Asbridge  et  al,  1968;  Scudder  et 
al,  1973;  Lin  et  al,  1974].  The  fact  that  the  enhanced  plasma  and  magnetic  turbulence 
and  their  associated  energetic  ions  are  similar  in  both  the  upstream  and  downstream 
regions  led  to  the  suggestion  of  the  same  bow  shock  source  region  [e.g..  West  and  Buck, 
1976',  Asbridge  et  al,  1978].  The  occurrence  rate  of  enhanced  energetic  ion  events 
increases  with  decreasing  Obu  both  upstream  aind  downstream  from  the  bow  shock  [e.g., 
Bonifazi  and  Moreno,  1981;  Mitchell  and  Roelof,  1983;  Crooker  et  al,  1981].  It  is  widely 
accepted  that  solar  wind  ions  undergo  first-order  Fermi  acceleration  by  scattering  back 
and  forth  across  the  quasi-parallel  bow  shock  in  the  turbulent  regions  upstream  and 
downstream  from  the  shock  [Lee,  1982;  Ellison,  1985].  In-situ  plasma  and  magnetic  field 
observations  have  confirmed  theoretical  predictions  of  the  coupling  between  the  particles 
and  waves  [Mobius  et  al,  1987;  Trattner  et  al,  1994].  Energy  spectra  of  energetic  ions 
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from  computer  simulations  are  also  in  a  very  good  agreement  with  observations  both 
upstream  and  downstream  from  the  quasi-parallel  shock  [Ellison  et  ai,  1990].  Bursts  of 
magnetospheric  ions  are  also  occasionally  observed  in  the  upstream  region  but  they  have 
a  much  harder  spectrum  than  those  of  bow  shock  accelerated  solar  wind  ions  [Mobius  et 
ai,  1986]. 

Ion  composition  measurements  have  played  a  critical  role  in  determining  the  relative 
contribution  of  bow  shock  accelerated  ions  and  leakage  of  magnetospheric  ions  to  the 
magnetosheath  energetic  ions  [e.g.,  Fuselier  et  al,  1991].  The  angular  distribution  of 
bow  shock  accelerated  ions  immediately  downstream  from  the  quasi-parallel  bow  shock 
is  anisotropic  with  the  maximum  flux  toward  the  magnetopause  [e.g.,  Scholer  et  ai, 
1989;  Ellison  et  al.,  1990].  They  contribute  to  the  majority  of  energetic  H'*'  (>  65%)  and 
almost  all  the  energetic  He"^^  in  the  magnetosheath  downstream  from  the  quasi-parallel 
shock  [Fuselier,  1994].  Leakage  of  magnetospheric  ions  plays  a  minor  role  in  this  region 
but  accounts  for  nearly  all  the  energetic  protons  in  the  plasma  depletion  layer  [Fuselier, 
1992],  a  thin  layer  adjacent  to  the  magnetopause  just  outside  the  magnetosphere  in  the 
magnetosheath  [Zwan  and  Wolf,  1976]. 

Recent  observations  of  cusp  energetic  particles  (CEP’s)  suggest  that  they  are  from 
the  bow  shock  accelerated  ions  [Chang  et  ai,  1998].  If  this  hypothesis  is  true,  one  would 
expect  the  appearance  of  bow  shock  accelerated  ions  in  the  magnetosheath  upstream 
from  the  cusp.  In  this  paper,  we  present  evidence  of  magnetosheath  energetic  ions 
upstream  from  the  cusp  accelerated  at  the  quasi-parallel  bow  shock  using  plasma  and 
magnetic  field  data  observed  by  the  NASA  GGS/POLAR  and  WIND  spacecraft  during 
a  magnetic  storm  event  on  May  4,  1998.  In  contrast  to  previous  statistical  work  [e.g., 
Crooker  et  ai,  1981;  Fuselier,  1994],  it  is  demonstrated  for  the  first  time  a  causal  relation 
between  the  magnetosheath  energetic  ions  and  bow  shock  magnetic  geometry.  This 
storm  provides  a  great  opportunity  to  understand  the  source  of  magnetosheath  energetic 
ions  and  to  check  the  hypothesis  of  bow  shock  source  of  CEP’s  by  suppling  large-scale 
temporal  variations  in  the  energetic  ion  fluxes  which  should  be  correlated  with  upstream 
parameters  at  the  source  region  with  a  proper  time  delay.  With  the  high-time  resolution 
data  from  POLAR  and  WIND,  based  on  the  aforementioned  principle  we  are  able  to 
determine  not  only  the  origin  of  the  observed  magnetosheath  energetic  ions  but  also  the 
upstream  parameters  that  control  the  intensity  of  these  ions. 

The  rest  of  the  paper  is  organized  as  followed.  In  section  2,  we  present  the 
instrumentation;  section  3,  POLAR  observations  of  plasma  composition  and  data 
interval  selection  for  cross-correlation  analysis;  section  4,  solar  wind  propagation  time; 
section  5,  magnetosheath  energetic  ions  and  their  correlation  with  IMF  cone  angle; 
section  6,  interpretation  of  the  analysis  results;  and  section  7  summary  and  conclusions. 
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2.  Instrumentation 

Energetic  ion  composition  measurements  during  this  storm  event  were  from  the 
Magnetospheric  Ion  Composition  Sensor  (MICS)  of  the  Charge  and  Mass  Magnetospheric 
Ion  Composition  Experiment  (CAMMICE)  onboard  POLAR.  A  similar  detector  was 
flown  on  the  CRRES  satellite  [Wilken  et  ai,  1992].  The  MICS  sensor  sampled  2-D 
angular  distributions  and  covered  the  energy  range  from  1  to  193.4  keV/e  for  H+  and 
1  to  100.1  keV/e  for  and  41.1  to  193.4  keV/e  for  It  took  about  3.3 

min  to  complete  a  full  energy  sweep.  The  H+  data  were  the  total  ion  measurements 
assuming  response  in  the  Double  Coincidence  Rate  (DCR)  channel  of  the  MICS 
detector.  The  detector  viewed  perpendicular  to  the  spacecraft  spin  axis  which  was 
mainly  in  the  dawn-dusk  direction  during  this  event.  Energetic  electron  observations 
were  from  the  Imaging  Electron  Sensor  (lES)  of  the  Comprehensive  Energetic  Particle 
and  Pitch  Angle  Distribution  (CEPPAD)  onboard  POLAR  [Blake  et  ai,  1995].  The 
sensor  sampled  3-D  distributions  covering  the  energy  range  from  12  to  400  keV.  The 
lower  energy  part  of  the  ion  spectrum  was  provided  by  the  Hydra  spectrometer  [Scudder 
et  ai,  1995]  which  sampled  3-D  distributions  of  electrons  and  ions  {assuming  H"'") 
covering  the  energy  range  from  about  10  eV  to  19  keV.  The  electron  measurements  from 
Hydra  as  well  as  the  magnetic  field  measurements  from  the  Magnetic  Field  Experiment 
(MFE)  magnetometer  [Russell  et  al,  1995]  were  used  to  identify  regions  of  plasma 
sheet,  magnetopause  boundary  layer,  magnetosheath,  and  solar  wind  along  the  POLAR 
trajectory.  Upstream  solar  wind  and  IMF  data  were  acquired  with  the  Magnetic  Field 
Investigation  (MFI)  magnetometer  [Lepping  et  ai,  1995]  and  the  Faraday  cup  of  the 
Solar  Wind  Experiment  (SWE)  [Ogilvie  et  ai,  1995]  both  onboard  the  WIND  spacecraft. 

3.  Overview  of  POLAR  Observations 

On  May  4,  1998,  0500-1200  UT,  the  POLAR  spacecraft  was  moving  from  —0.3°  to 
74.1°  magnetic  latitude  and  from  4.4  to  9.0  Re  geocentric  distance  near  the  magnetic 
local  noon.  The  WIND  spacecraft  was  about  214  Re  upstream  from  the  Earth,  just 
north  of  the  Sun-Earth  line.  A  very  strong  interplanetary  shock  arrived  at  WIND  at 
about  0230  UT.  WIND  observed  large  variations  in  the  solar  wind  density  (~4-65  cm“^) 
and  IMF  strength  (~5-40  nT)  and  high  solar  wind  bulk  speed  ~800  km/s  for  the  next  9 
hours.  Because  of  high  solar  wind  dynamic  pressure  (as  large  as  65  nP),  magnetosphere 
was  severely  compressed  and  bow  shock  even  reached  the  POLAR  altitude  on  several 
occasions.  As  a  result  of  the  rapid  boundary  motion  due  to  the  dynamics  of  the  storm, 
POLAR  frequently  crossed  the  magnetopause  current  layer.  However,  POLAR  was  in 
the  high-latitude  dayside  magnetosheath  most  of  the  time  during  the  second  half  period 
of  the  storm  event.  (Note  only  relevant  POLAR  data  are  presented  here  to  address 
the  plasma  composition  and  magnetosheath  interval  selection  for  the  cross-correlation 
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analysis.  For  detailed  analysis  of  POLAR  observations  for  this  event  see  Scudder  et  al. 
(AND/OR  OTHERS?)  [1999d].) 

Plate  1  presents  from  top  to  bottom  panels  MFE  magnetic  field  magnitude, 
energy-time  spectrograms  for  CAMMICE/MICS  DCR-H"^,  Hydra 

ions,  CEPPAD/IES  electrons,  and  Hydra  electrons  from  0500  to  1200  UT.  The 
total  energy  coverages  for  H'*'  and  electron  are  ~0. 02-200  keV  and  ~0. 01-400  keV, 
respectively.  Before  0525  UT,  POLAR  was  in  the  plasma  sheet  characterized  by  a  large 
magnetic  field  (>  220  nT)  dominated  by  the  positive  R.  component,  and  relatively 
intense  energetic  (>  40  keV/e),  H+  (>  40  keV),  and  electrons  (>  20  keV)  but 

relatively  weak  lower-energy  (below  40  keV/e)  He"^^,  and  H"*'  fluxes.  POLAR 

then  moved  into  the  low-latitude  boundary  layer  (LLBL)  and  observed  similar  features 
before  it  exited  the  magnetopause  into  the  magnetosheath  at  about  0541  UT  [Russell 
et  ai,  1999].  Although  the  plasma  sheet  and  LLBL  are  both  on  closed  magnetic  field 
lines  with  both  feet  anchored  to  the  Earth,  electron  pitch  angle  distributions  there  are 
somewhat  different,  trapped  population  in  the  former  region  and  counter-streaming 
distribution  in  the  latter. 

By  contrast,  in  the  magnetosheath  from  ~0541  to  ~0610  UT  POLAR  observed 
disturbed  weak  magnetic  field,  enhanced  0>+“,  He^^,  H+  and  1-keV  ions  and  100-eV 
electrons  but  weak  and  energetic  electron  flux  (both  nearly  close  to  the  background 

values).  Then  POLAR  moved  into  the  LLBL  again.  During  the  early  part  of  the 
storm  event  from  0500  to  0730  UT,  POLAR  was  moving  back  and  forth  between  the 
magnetosheath  and  the  plasma  sheet/LLBL  across  the  magnetopause.  After  0730  UT, 
POLAR  was  in  the  magnetosheath  and  encountered  the  bow  shock  for  the  first  time 
at  about  0735  UT.  Due  to  the  arrival  of  a  brief,  very  strong  pressure  pulse  (solar  wind 
dynamic  pressure  ~65  nP)  bow  shock  moved  to  7.3  Re  at  38°  magnetic  latitude  and 
11  o’clock  magnetic  local  time.  POLAR  briefly  stayed  upstream  from  the  shock  and 
observed  very  cold  solar  wind  ions  and  electrons  as  shown  in  the  Hydra  ion  and  electron 
spectrograms  in  Plate  1.  It  encountered  bow  shock  several  times  after  that. 

The  bow  shock  crossings  are  more  obviously  illustrated  in  Figure  1  which  shows  the 
magnitude  of  the  magnetic  field  from  MFE  and  electron  density  from  Hydra  covering 
0730-1000  UT.  Bow  shock  crossings  are  identified  by  the  simultaneous,  abrupt  jumps 
in  the  two  quantities.  In  addition  electron  bulk  speed  and  temperature  jumps  also 
occurred  corresponding  to  the  magnetic  field  and  electron  density  jumps.  As  POLAR 
moved  from  the  magnetosheath  across  the  strong  shock  into  the  upstream  region,  a 
sharp  increase  in  the  electron  bulk  speed  and  a  sharp  decrease  in  the  magnetic  field, 
electron  temperature  and  density  occurred,  and  vice  versa  [Scudder  et  al.  (AND/OR 
OTHERS?),  1999c].  Several  upstream  regions  from  the  bow  shock  are  identified  and 
marked  as  the  shaded  regions  in  Figure  1.  As  the  solar  wind  dynamic  pressure  quickly 
reduced,  the  bow  shock  and  magnetopause  retreated  toward  their  normal  positions  and 
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Plate  1.  Magnetic  field  and  plasma  measurements  from  four  instruments  onboard  PO¬ 
LAR  from  5  to  12  UT  on  May  4,  1998:  from  top  to  bottom,  MFE  magnetic  field  mag¬ 
nitude,  energy-time  spectrograms  for  CAMMICE/MICS  0<+^  0>+2^  He+2,  DCR-H+, 
Hydra  ions,  CEPPAD/IES  electrons,  and  Hydra  electrons.  Distinct  features  were  appar¬ 
ent  as  POLAR  traversed  various  regions  between  the  plasma  sheet  and  upstream  region 
of  the  bow  shock. 
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POLAR  crossed  the  magnetopause  into  the  magnetopause  boundary  layer  at  about  0755 
UT.  Electron  temperature  then  increased  as  expected. 

The  best  signature  however  for  distinguishing  the  boundary  layer  and  magnetosheath 
region  is  the  electron  distribution  since  two  regions  have  distinct  magnetic  topology, 
one  for  a  closed  magnetic  field  line  and  the  other  for  a  magnetosheath  field  line  which 
is  connected  to  IMF  and  may  or  may  not  be  connected  to  a  geomagnetic  field  line. 

As  an  example  shown  in  Figure  2a,  electrons  commonly  show  a  counter-streaming 
distribution  in  the  boundary  layer,  indicating  a  closed  magnetic  topology.  Another 
very  useful  diagnostic  of  the  magnetopause  current  layer  as  a  rotational  discontinuity  is 
the  generalized  Walen  test  [Scudder  et  ai,  1999a].  Magnetopause  crossings  identified 
using  this  technique  are  consistent  with  the  findings  based  on  the  electron  distributions 
[Scudder  et  ai  (REFERENCE?),  1999b].  After  0755  UT,  POLAR  frequently  crossed  the 
magnetopause  before  exiting  the  boundary  layer  at  ~0836  UT.  After  that,  except  several 
brief  upstream  regions  from  the  bow  shock,  POLAR  was  in  the  magnetosheath  until 
the  end  of  the  magnetic  storm  at  about  1200  UT.  The  common  electron  distribution  in 
the  magnetosheath  is  nearly  isotropic,  accompanied  by  a  weak  field-aligned  beam  at  low 
energy  as  shown  in  Figure  2b. 

As  shown  in  Plate  1,  ion  composition  of  energetic  ions  (O^"*"^,  O^'''^,  He"*"^,  and 
H"*')  shows  a  distinct  relative  abundance  in  the  plasma  sheet  and  magnetosheath. 
Magnetosheath  energetic  ions  are  likely  of  the  solar  wind  origin.  During  this  storm 
event,  the  magnetosheath  energetic  ion  fluxes  show  lots  of  variations  (see  for  example 
50  keV/e  ions  in  the  third  to  the  fifth  panel  of  Plate  1).  Energetic  ion  fluxes  observed 
upstream  in  the  solar  wind  by  the  3D  Plasma  instrument  [Lin  et  ai,  1995]  onboard 
WIND  on  the  other  hand  were  quite  uniform.  In  addition,  these  fluxes  are  about  1 
to  3  orders  of  magnitude  lower  than  the  magnetosheath  energetic  ion  fluxes  observed 
by  POLAR.  Therefore,  the  solar  energetic  particles  can  not  be  a  direct  source  of 
the  POLAR  observations.  To  understand  the  cause  of  the  variations  and  ultimately 
the  origin  of  these  magnetosheath  energetic  ions,  we  choose  the  later,  long  period  of 
magnetosheath  interval  (0840-1200  UT  excluding  upstream  intervals)  for  systematic 
analysis. 

4.  Solar  Wind  Propagation  Time 

If  the  source  of  the  magnetosheath  energetic  ions  is  at  a  distance,  the  timing  for 
the  transport  is  a  crucial  factor.  Since  their  composition  suggests  a  solar  wind  source, 
solar  wind  propagation  time  becomes  the  next  obvious  quantity  to  examine.  The 
propagation  time  can  be  estimated  by  four  different  methods  utilizing  the  plasma  and 
field  measurements  from  WIND  and  POLAR.  The  first  simple  estimate  is  to  use  the 
measured  solar  wind  bulk  speed  from  WIND/SWE  which  was  quite  uniform  ~770  km/s 
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Figure  1.  POLAR/MFE  magnetic  field  magnitude  and  Hydra  electron  density  from 
0730  to  1000  UT.  POLAR  was  mainly  in  the  dayside  magnetosheath  and  briefly  crossed 
the  bow  shock  into  the  upstream  region  (shaded)  several  times  as  indicated  by  the  si¬ 
multaneous,  abrupt  jumps  in  the  two  quantities. 
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(a)  05-04-1998  07:59:31,73-07:59:44.23 


Figure  2.  Contours  of  constant  phase  space  density  for  electrons  from  POLAR/Hydra  in 
the  plane  of  velocity  components  perpendicular  and  parallel  to  the  magnetic  field,  showing 
(a)  a  counter-streaming  distribution  at  energies  below  ~500  eV  in  the  magnetopause 
boundary  layer  and  (b)  a  nearly  isotropic  distribution  with  a  weak  field-aligned  beam  at 
~60  eV  in  the  magnetosheath. 
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from  0650  to  0730  UT.  Estimated  propagation  time  is  28  min  for  solar  wind  traveling  a 
distance  of  204  Rg  from  WIND  to  bow  shock. 

The  POLAR  magnetic  field  profile  observed  in  the  upstream  region  of  the  bow 
shock  from  0735:47  to  0738:03  UT  is  unique  to  the  WIND  IMF  data  set  so  that  a 
comparison  of  IMF  measurements  from  two  spacecraft  can  give  a  fairly  accurate  estimate 
of  the  solar  wind  propagation  time.  Figure  3  presents  the  magnetic  field  measurements 
from  POLAR/MFE  and  WIND/MFI  offset  27.3  min.  The  scales  of  all  the  vertical  axes 
are  set  to  facilitate  the  comparison  and  therefore  most  of  the  POLAR  measurements  are 
out  of  scale.  In  the  upstream  region,  and  components  from  POLAR  both  turned 
from  positive  to  negative  at  about  the  same  time.  An  identical  feature  was  also  observed 
by  WIND  about  27.3  min  ahead.  All  the  three  components  and  magnitude  of  IMF  from 
POLAR  and  WIND  (shaded  region)  match  very  well,  with  the  exception  of  the  By 
components  which  show  different  trends  but  maintain  the  same  sign.  Within  2  hours 
ahead  of  this  POLAR  upstream  interval,  only  one  match  for  all  three  IMF  components 
occurs  in  the  WIND  data  set.  Other  intervals  have  at  best  a  match  in  one  component 
only.  The  solar  wind  propagation  time  from  WIND  to  POLAR  just  upstream  from  the 
bow  shock  is  thus  27.3  min  consistent  with  the  estimate  from  the  solar  wind  bulk  speed. 
For  the  time  interval  of  interest,  0810-1130  UT,  corresponding  to  the  interval  of  the 
POLAR  magnetosheath  energetic  ion  observations,  the  solar  wind  bulk  speed  decreased 
to  about  755  km/s  which  gives  an  estimate  of  the  propagation  time  of  ~29  min  from 
WIND  to  bow  shock.  Additional  several  minutes  are  required  for  plasmas  to  propagate 
from  the  shock  to  POLAR  in  the  magnetosheath  because  bow  shock  decelerates  the 
flow. 

Since  magnetosheath  magnetic  field  lines  are  connected  to  the  IMF ,  magnetosheath 
field  and  IMF  to  some  extent  would  be  correlated  despite  magnetic  field  being  modified 
by  the  electric  current  at  bow  shock.  Figure  4  presents  the  three  magnetic  field 
components  observed  by  POLAR  from  0840  to  1200  UT  and  by  WIND  with  a  lag 
indicated  in  each  panel.  These  lags  are  determined  by  a  cross-correlation  analysis  of 
the  POLAR  and  WIND  observations  during  the  above  interval  excluding  the  upstream 
regions  from  bow  shock  (shaded  regions  in  the  figure).  Correlation  coefficient  for  each 
magnetic  field  component  is  calculated  for  a  lag  from  0  to  60  min.  All  three  correlation 
coefficient  curves  show  one  unique  peak  with  the  correlation  coefficients  (0.94,  0.46, 
0.73)  and  time  lags  (35,  34,  37)  min  for  B^,  By,  and  B,,  respectively.  Strong  correlation 
appears  in  B^,  good  correlation  in  Bz  but  weak  correlation  in  By,  as  illustrated  in 
the  figure.  From  these  results,  the  estimated  plasma  propagation  time  from  WIND  to 
POLAR  is  approximately  35-37  min  which  is  consistent  with  the  estimate  using  the 
solar  wind  bulk  speed  and  magnetosheath  plasma  flow  speed. 

Another  estimate  of  the  propagation  time  is  given  by  the  cross-correlation  analysis 
of  the  solar  wind  density  and  magnetosheath  plasma  density.  Figure  5  presents  electron 
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Figure  3.  Magnitude  and  three  GSM  components  of  magnetic  field  observed  by  PO- 
LAR/MFE  (thin)  from  0700  to  0750  UT  and  WIND/MFI  (thick)  in  the  solar  wind  offset 
27.3  min.  The  IMF  profile  from  POLAR  just  upstream  from  the  bow  shock  (shaded 
interval)  matches  the  WIND  observations  very  well.  The  scales  of  the  vertical  axes  are 
set  for  an  easy  comparison  of  IMF  profiles.  Therefore,  MFE  measurements  are  off  scale 
most  of  the  time. 
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Figure  4.  GSM  magnetic  field  components  in  the  magnetosheath  from  POLAR/MFE 
(thick)  from  0840  to  1200  UT  and  IMF  components  from  WIND/MFI  (thin)  offset  by  a 
At.  The  At  in  each  panel  is  the  time  delay  associated  with  the  peak  correlation  between 
the  magnetosheath  magnetic  field  and  IMF.  Intervals  of  shaded  areas  are  excluded  in  the 
cross-correlation  analysis.  Very  strong  correlation  appears  in  5^,.,  good  correlation  in 
but  weak  correlation  in  By.  Estimated  solar  wind  propagation  time  is  ~  35-37  min. 
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density  from  POLAR/Hydra  from  0840  to  1200  UT  and  solar  wind  proton  density  from 
WIND/SWE  offset  33  min.  The  upstream  regions  (shaded)  are  again  excluded  in  the 
analysis.  As  expected,  magnetosheath  plasma  density  is  strongly  correlated  with  the 
solar  wind  density  with  a  proper  time  delay.  The  correlation  coefficient  as  a  function 
of  the  lag  from  0  to  60  min  once  again  demonstrates  a  unique  peak  value  of  0.85  at  a 
lag  of  33  min,  yielding  another  estimate  of  the  solar  wind  propagation  time.  This  peak 
however  is  not  as  sharp  as  those  in  the  analysis  of  the  magnetosheath  magnetic  field 
and  IMF  and  5,,..  components  presented  above.  Therefore,  there  is  a  somewhat  larger 
uncertainty  (±3  min)  in  the  estimate.  Nevertheless,  the  estimated  propagation  time  is 
consistent  with  the  one  based  on  the  magnetic  field  comparison. 

5.  Magnetosheath  Energetic  Ions 

Figure  6  depicts  energy  spectra  of  H''",  He"^',  and  from  Hydra  and  CAMMICE 
in  the  magnetosheath,  averaged  over  an  interval  from  ~0840  to  ~1200  UT  excluding 
the  solar  wind  intervals.  Ion  fluxes  are  plotted  against  the  ion  energy  per  charge.  In 
this  format,  spectral  shapes  for  all  three  species  are  similar.  In  particular,  they  are 
nearly  identical  for  H"*"  and  with  energy  above  20  keV /e.  Some  differences  in 
are  due  to  very  low  count  rates  which  create  much  larger  errors  in  the  counting 
statistics.  It  gets  worse  as  energy  gets  higher.  Nevertheless,  all  three  spectra  are 

continuous  with  a  spectral  break  at  about  41.1  keV/e.  These  spectra  resemble  those 
upstream  and  downstream  from  the  cpiasi-parallel  bow  shock  (cf.  Figures  3  and  4  of 
Mohius  et  al.  [1987]  and  Figure  2  of  Ellison  et  al.  [1990]).  It  is  noted  that  there  are 
substantial  variations  in  the  ion  fluxes  during  this  magnetosheath  interval.  Ion  spectra 
of  the  highest  time  resolution  (3.3  min)  would  move  up  and  down  along  the  vertical  axis 
in  Figure  6  as  time  steps  through  each  interval  of  full  energy  sweep  within  the  above 
magnetosheath  interval  but  their  spectral  shapes  remain  similar. 

To  examine  the  hypothesis  of  the  observed  magnetosheath  energetic  ions,  especially 
ions  with  energy  above  the  spectral  break,  originated  from  the  bow  shock  accelerated 
ions,  a  cross-correlation  analysis  of  the  magnetosheath  H"*"  flux  integrated  from  41.1  to 
193.4  keV/e  and  the  IMF  cone  angle  6bx  (the  angle  between  the  IMF  and  Sun-Earth 
line)  was  performed.  Figure  7a  presents  the  integrated  H'*'  flux  at  the  highest  time 
resolution  from  0840  to  1200  UT  and  9bx  offset  36  min.  Both  quantities  show  large 
temporal  variations,  the  maximum  to  minimum  flux  ratio  ~100  and  6bx  ~  10°-90°.  In 
addition,  H'*'  flux  increases  as  9bx  decreases,  and  vice  versa.  Two  c|uantities  are  strongly 
anti-correlated  with  a  correlation  coefficient  of  —0.82. 

The  correlation  coefficient  for  the  energetic  ion  flux  and  6bx  during  this 
magnetosheath  interval  in  general  depends  on  the  time  lag  for  6bx^  shown  in  Figure  7b. 
As  the  time  lag  increases  from  0  to  60  min,  correlation  coefficient  starts  at  about 
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UT 

Figure  5.  POLAR/Hydra  electron  density  in  the  magnetosheath  (thick)  from  0840  to 
1200  UT  and  the  solar  wind  proton  density  from  WIND/SWE  (thin)  offset  33  min.  Two 
quantities  are  strongly  correlated.  The  33-min  lag  associated  with  the  best  correlation 
gives  an  estimate  of  the  solar  wind  propagation  time.  Intervals  of  shaded  areas  are 
excluded  in  the  cross^correlation  analysis. 
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Figure  6.  Energy  spectra  of  magnetosheath  H"*",  He"^",  and  from  POLAR/Hyclra 

and  CAMMICE  averaged  over  the  interval  ~0840-1200  UT  excluding  the  solar  wind 
intervals.  Spectral  shapes  are  similar,  ordered  by  the  ion  energy  per  charge.  They  are 
nearly  identical  for  H'*'  and  He'*'-  with  energy  above  20  keV/e.  Some  differences  in 
are  due  to  low  count  rate.  Spectra  are  continuous  with  a  spectral  break  at  ~41.1  keV/e 
(vertical  dashed  line). 
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0,  monotonically  decreases,  reaches  the  peak  value  —0.82  at  36-min  lag,  and  then 
monotonically  increases  toward  0.  This  36-min  lag  associated  with  the  best  correlation 
is  consistent  with  the  estimated  solar  wind  propagation  time  and  other  time  estimates 
(see  discussions  below). 

Although  all  three  ion  spectra  in  Figure  6  are  continuous  with  a  spectral  break 
at  the  same  energy,  lower-  and  higher-energy  ions  demarcated  by  this  energy  have 
distinct  properties  with  respect  to  6 bx-  Similar  to  the  format  in  Figure  7b,  correlation 
coefficients  for  the  magnetosheath  ion  flux  integrated  from  Ei  to  E2  and  6bx  plotted 
against  the  time  lag  from  0  to  60  min  are  presented  in  Plate  2  for  three  ion  species  (H"*", 
He'*'^,  where  Ei  steps  through  the  detector  energy  channel  from  1.0  to  100.1 

keV/e  and  E2  is  fixed  at  the  highest  detector  energy  channel,  193.4  keV/e  for  H'*'  and 
He+^  and  100.1  keV/e  for  These  correlation  curves  are  color  coded  according  to 

their  EiS  given  in  the  color  bar.  The  first  most  remarkable  result  is  that  all  the  curves 
at  F'l  >  41.1  keV/e  (in  blue)  from  panels  a  and  b  are  nearly  identical.  In  addition, 
these  curves  have  a  unique  peak  correlation  coefficient  ~  —0.8  at  a  time  lag  of  ~36-37 
min,  suggesting  that  magnetosheath  energetic  (>  41.1  keV/e)  H'''  and  He+2  ion  flux  are 
strongly  anti-correlated  with  6bx- 

It  is  even  more  striking  that  correlation  curves  for  H"''  and  He"^^  begin  to  change 
drastically  by  including  ion  fluxes  from  one  or  two  energy  channels  below  41.1  keV/e. 
They  become  much  flatter  and  closer  to  the  line  of  null  correlation  as  shown  by  the  green 
curves  in  panels  a  and  b.  This  result  suggests  a  critical  energy  for  the  anti-correlation 
between  the  magnetosheath  energetic  ion  flux  and  9bx-  As  Ei  continues  to  decrease 
toward  1  keV /e,  the  correlation  curves  for  H+  and  He^^  begin  to  show  some  positive 
correlation  and  then  move  toward  weak  correlation  as  indicated  by  the  yellow  and  red 
curves.  Because  energetic  ions  have  very  low  count  rates,  correlation  curves  are 

rather  irregular  as  shown  in  panel  c.  Nevertheless,  these  curves  show  a  trend  similar  to 
those  for  H+  and  He+^.  Furthermore,  at  lower  energies  ions  gain  enough  counting 

statistics  so  that  correlation  coefficient  curves  have  more  significance.  These  curves  are 
nearly  identical  to  those  of  He+^  at  the  same  energies  (cf.  red  curves  in  panels  b  and  c). 

To  illustrate  the  strong  energy  dependence  in  the  correlation  coefficient,  the  peak 
correlation  coefficient  is  plotted  against  Ei  for  all  three  ion  species  as  shown  in 
Figure  8a.  It  is  clearly  demonstrated  that  a  critical  energy  Ec  has  to  be  reached  for 
ion  flux  to  be  anti-correlated  with  9bx-  This  critical  energy  equals  to  41.1  keV/e  for  all 
three  ion  species  and  is  exactly  identical  to  the  one  at  the  spectral  break  suggested  by 
Figure  6.  Furthermore,  r^’s  are  nearly  constant  for  Ei  above  E^  ~  —0.82  for  both 
and  He"*"^  and  rvj  —  0.71  for  0^+^.  However  when  Ei  is  above  70  keV/e,  begins  to 
move  toward  0  due  to  less  significance  in  the  ion  counting  statistics.  This  trend  is  more 
obvious  for  higher-energy  since  their  fluxes  are  nearly  close  to  the  background 

value.  Nonetheless,  the  sharp  transition  from  strong  anti-correlation  to  moderate 
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(b)  Correlation  Coefficient  vs  Time  Lag 


Figure  7.  (a)  Energetic  (41.1  <  E  <  193.4  keV)  H"*"  ion  flux  (thick  line)  from  PO- 
LAR/CAMMICE  in  the  interval  0840-1200  UT  and  ^bx  (thin  line)  from  WIND/MFI 
offset  36  min.  Both  show  large  temporal  variations  and  are  strongly  anti-correlated,  (b) 
The  correlation  coefficient  is  a  function  of  the  time  lag.  As  time  lag  increases  from  0  to 
60  min,  the  correlation  coefficient  starts  at  about  0,  monotonically  decreases,  reaches  the 
peak  value  —0.82  at  a  lag  of  36  min,  and  then  monotonically  increases  toward  0.  Note 
that  the  shaded  intervals  are  excluded  in  the  calculation  of  the  correlation  coefficient. 
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Plate  2.  Correlation  coefficient  curves  of  the  magnetosheath  ion  flux  {Ei  <  E  <  E2) 
and  6bx  for  three  ion  species,  (a)  H+,  (b)  He'''^,  and  (c)  Curves  are  color  coded 

according  to  their  Ei's  given  in  the  color  bar  with  E2  equal  to  193.4  keV/e  for  and 
He"*"^  and  101.1  keV/e  for  All  the  curves  of  H''"  and  He'*"^  with  Ei  >  41.1  keV/e  (in 

blue)  are  nearly  identical,  showing  a  strong  anti-correlation  between  the  ion  flux  and  Qbx 
offset  36-37  min  (peak  correlation  coefficient  ~  —0.8).  Below  41.1  keV/e,  curves  sharply 
shift  toward  weakly  positive  correlation  (in  red).  Similar  trends  also  appear  in 
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positive  correlation  just  below  Ef.  is  vividly  illustrated. 

The  time  lag  Ato  associated  with  the  peak  correlation  can  be  quite  different  for 
different  energies  as  shown  in  Figure  8b.  However  when  r„  is  very  small  {<  —0.7), 

At,,  is  nearly  the  same  for  all  three  ions  at  about  36-37  min.  This  At,,  suggests  a 
unic|ue  signal  transmission  time  independent  of  ion  energy  and  species.  Results  of  the 
cross-correlation  analysis  for  all  three  ions  have  demonstrated  that  magnetosheath 
energetic  ions  with  energy  above  the  spectral  break  are  strongly  anti-correlated  with 
offset  by  a  lag  consistent  with  the  solar  wind  propagation  time  and  other  time  estimates 
that  will  be  discussed  later.  This  correlation  cannot  be  incidental.  Results  of  r„  and  At„ 
for  lower-energy  [Ei  <  10  keV/e)  ions  shown  in  Figure  8  on  the  other  hand  can  only 
indicate  average  results  for  lower-energy  ions  and  will  be  discussed  more  later. 

Since  the  observed  magnetosheath  energetic  ion  fluxes  are  controlled  by  Ob-,-.,  one 
would  expect  they  flow  antisunward  away  from  the  bow  shock.  In  this  storm  event, 
the  magnetosheath  magnetic  field  was  most  of  the  time  close  to  the  POLAR’s  spin 
axis  so  that  CAMMICE  detector  only  covered  a  limited  range  of  pitch  angle  around 
90°.  However  during  the  interval  from  ~0850  to  1000  UT,  magnetosheath  energetic 
ions  have  a  nearly  full  coverage  of  pitch  angle  that  can  deduce  their  flow  direction. 
Figure  9  presents  the  normalized  angular  distribution  of  energetic  ion  flux  at  the 
highest  time  resolution  for  H"*",  He"^^  (both  41.1  <  E  <  193.4  keV/e),  and 
(41.1  <  E  <  101.1  keV/e)  in  the  spacecraft  frame  for  three  different  intervals  in  the 
magnetosheath.  Because  the  plasma  flow  speed  in  the  POLAR  spin  plane  is  about  300 
km/s  corresponding  to  ~0.5  keV/e  in  energy  for  H"^  and  He'*'^,  correction  by  the  plasma 
flow  for  the  above  energetic  ion  distributions  is  expected  to  be  very  small.  Therefore, 
Figure  9  is  a  very  good  approximation  of  the  energetic  ion  distributions  in  the  rest  frame 
of  plasma. 

As  shown  in  Figure  9b  and  c,  almost  all  the  H"''  and  He'^'^  in  the  magnetosheath 
proper  and  near  the  bow  shock  have  pitch  angles  less  than  90°  with  the  maximum  flux 
occurring  at  about  45°  pitch  angle.  Therefore  these  ions  show  a  strong  flow  along  the 
magnetic  field  direction.  During  the  interval  covered  in  Figure  9a,  POLAR  was  in  the 
magnetosheath  close  to  the  magnetopause  ajid  energetic  ions  with  pitch  angles  within 
30°  to  the  magnetic  field  line  were  not  sampled.  Nevertheless  the  observed  partial 
angular  distributions  for  H"''  and  He"^^  in  this  region  are  consistent  with  those  observed 
in  the  magnetosheath  proper  and  near  the  bow  shock.  As  mentioned  several  times 
before,  energetic  ions  are  tenuous  that  is  more  apparent  in  this  diagram.  However, 
their  angular  distributions  show  a  trend  similar  to  those  of  the  other  two  species. 

Within  the  above  magnetosheath  interval  (0850-1000  UT),  the  angular  distributions 
of  energetic  ions  are  basically  the  same  as  those  presented  in  Figure  9,  showing  a 
predominant  flow  along  the  magnetic  field  direction.  Meanwhile  the  magnetosheath 
magnetic  field  component  turned  negative  and  turned  positive  after  0840  UT 
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Figure  8.  (a)  Peak  correlation  coefficient  as  a  function  of  energy  Ei  for  three  ion 
species.  For  Ei  >  Ec  =  41.1  keV/e,  r^’s  are  nearly  the  same,  ~  —0.82  for  both  H+ 
and  He'*'^  and  ~  —0.71  for  0^+^.  For  Ei  above  70  keV/e,  r„’s  turn  toward  0  due  to 
lower  count  rate  especially  in  A  sharp  transition  from  strong  anti-correlation  to 

weakly  positive  correlation  occurs  at  energy  immediately  below  E^  This  critical  energy 
Ec  happens  to  be  the  energy  at  the  spectral  break  suggested  by  Figure  6.  (b)  Time  lag 
associated  with  the  peak  correlation  as  a  function  of  energy  Ei.  In  general  Ato  can 
be  quite  different.  All  the  points  of  large  negative  r^’s  (<  -0.7)  with  Ei  >  Ec  however 
have  Ato  =  36-37  min. 
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and  maintained  the  same  polarity  until  ~1002  UT  as  shown  in  Figure  4.  Therefore 
the  observed  magnetosheath  energetic  ions  during  the  above  interval  were  streaming 
antisunward  and  northward  away  from  the  bow  shock. 

6.  Discussions 

In  the  previous  section  we  have  presented  plasma  and  magnetic  field  measurements 
from  the  NASA  GGS/POLAR  and  WIND  satellites  for  the  May  4,  1998  storm  event. 
Here  we  discuss  origin  of  the  observed  magnetosheath  energetic  ions  based  on  possible 
explanations  in  the  literature.  In  general  two  major  sources  are  the  solar  wind  and 
magnetospheric  plasmas  but  the  acceleration  region  and  the  transport  can  be  quite 
different  for  each  of  the  sources.  Energetic  magnetospheric  ions  in  the  plasma  sheet 
on  closed  field  lines  can  leak  into  the  magnetosheath  [e.g.,  Sibeck  et  ai,  1987].  Their 
intensity  is  more  related  to  substorm  activities  in  the  magnetotail.  Energetic  ions 
previously  on  closed  field  lines  can  escape  the  magnetosphere  into  the  magnetosheath 
along  open  field  lines  either  by  gradient  and  curvature  drift  across  field  lines  or 
magnetic  reconnection  [e.g.,  Scholer  et  ai,  1981].  Solar  energetic  particles  can  also 
directly  contribute  to  the  magnetosheath  energetic  ions.  Diffuse  ions  accelerated  at  the 
quasi-parallel  bow  shock  can  be  another  source  [e.g.,  West  and  Buck,  1976;  Fuselier  et 
ai,  1991]. 

The  composition  of  the  observed  magnetosheath  energetic  ions  suggest  they  are 
ultimately  from  the  solar  wind  not  from  the  ionosphere.  It  is  a  question  whether  keV 
ions  of  solar  wind  origin  entered  the  magnetosphere,  convected  to  the  magnetotail,  were 
accelerated  there  and  injected  to  the  dayside,  and  then  escaped  the  magnetosphere. 
Because  the  dayside  plasma  sheet  and  magnetosheath  traversed  by  POLAR  during 
this  event  show  distinct  energetic  ion  composition,  leakage  of  magnetospheric  ions  is 
not  likely  a  cause  of  the  magnetosheath  energetic  ions.  In  addition,  POLAR  was  on 
magnetosheath  field  lines  which  were  most  of  the  time  connected  to  the  solar  wind,  not 
the  open  magnetosphere  (cusp,  mantle,  lobe,  etc.).  Escape  of  energetic  ions  through  the 
open  field  lines  is  not  a  possible  explanation  for  the  observations.  Most  importantly 
magnetospheric  source  would  have  a  difficulty  to  explain  the  fact  that  magnetosheath 
energetic  ion  flux  of  solar  wind  origin  are  anti-correlated  with  the  IMF  cone  angle  {6bx)- 

Solar  energetic  ions  observed  by  WIND  are  much  less  intense  than  the 
magnetosheath  energetic  ions  during  this  event.  Results  of  cross-correlation  analysis 
reveal  that  both  ion  fluxes  are  slightly  anti-correlated.  Therefore  solar  energetic  ions  are 
not  a  direct  source  of  magnetosheath  energetic  ions.  However,  an  acceleration  region 
at  the  bow  shock  can  account  for  ions  accelerated  to  hundreds  of  keV.  As  noted  a 
long  time  ago,  energetic  ion  spectra  downstream  from  the  shock  are  nearly  identical 
to  those  in  the  upstream  wave  region  (see  for  example  Figure  10  of  West  and  Buck 
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Figure  9.  Normalized  angular  distribution  of  energetic  ion  flux  for  H'*',  He'*'^  (both 
41.1  <  E  <  193.4  keV/e),  and  0^"^^  (41.1  <  E  <  101.1  keV/e)  in  the  spacecraft  frame 
for  three  intervals  in  the  magnetosheath,  (a)  near  the  magnetopause  boundary  layer, 
(b)  in  the  magnetosheath  proper,  and  (c)  near  the  bow  shock.  and  ions  are 
streaming  along  the  magnetic  field  direction.  Note  that  pitch  angles  within  30°  to  the 
magnetic  field  line  in  (a)  are  not  sampled.  Tenuous  ions  are  apparent  in  this 

diagram  but  they  show  a  trend  consistent  with  the  other  two  species. 
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[1976]).  Ions  undergo  first-order  Fermi  acceleration  by  scattering  back  and  forth  across 
the  quasi-parallel  bow  shock  in  the  turbulent  regions  upstream  and  downstream  from 
the  shock.  This  mechanism  is  most  efficient  for  highly  charged  ions  [Lee,  1982].  The 
magnetosheath  energetic  ion  spectra  observed  by  POLAR  have  shown  such  a  property. 

Results  of  the  cross-correlation  analysis  for  the  magnetosheath  ion  flux  and  Obt.  are 
very  striking.  In  order  to  lower  the  statistical  error  especially  for  higher-energy  ions,  we 
have  done  this  analysis  using  the  integrated  ion  number  flux  instead  of  the  differential 
number  flux  at  each  MICS  energy  channel.  Because  of  the  steep  drop  in  the  energetic 
tail  of  the  ion  spectra,  ion  flux  detected  at  the  lowest  energy  channel  in  the  integration 
energy  range  contributes  most  to  the  total  flux.  Therefore  the  integrated  flux  for  ion 
energy  above  ~10  keV/e  to  some  extent  truly  reflects  the  behavior  of  the  differential 
flux  at  the  lowest  energy  as  demonstrated  by  the  sharp  transition  in  Figure  8.  In  fact, 
the  correlation  coefficient  curves  for  ion  differential  number  flux  at  energy  above  ~10 
keV/e  are  very  similar  to  the  corresponding  ones  for  the  integrated  flux  presented  in 
Plate  2.  On  the  other  hand,  for  ion  energy  below  10  keV/e  the  ion  differential  energy 
flux  at  the  lowest  energy  does  not  dominate  the  total  flux  because  the  ion  flux  at 
this  energy  is  not  large  enough  to  compensate  the  reduced  energy  in  the  integration. 
Thus,  results  of  correlation  analysis  including  these  ions  presented  here  only  represent 
average  results  rather  than  results  for  ions  detected  at  the  lowest  energy  channel  within 
the  integration.  The  correlation  curves  for  the  differential  ion  flux  at  energy  below  10 
keV/e  become  more  erratic  showing  a  trend  toward  null  correlation  with  decreasing 
energy,  unlike  the  smooth,  asymptotic  behavior  shown  in  Figure  8.  This  result  of  little 
correlation  is  expected  because  keV  magnetosheath  ions  mainly  come  from  the  subsolar 
bow  shock  along  streamlines  and  the  intensity  of  solar  wind  ions  there  is  not  related  to 
the  IMF  or  bow  shock  geometry. 

The  most  remarkable  feature  in  the  results  of  the  cross-correlation  analysis  is  that 
the  correlation  coefficient  is  energy  dependent.  Correlation  coefficients  as  a  function  of 
ion  energy  show  a  sharp  transition  at  the  same  energy  (precisely  energy  per  charge)  for 
all  three  ion  species.  Above  this  critical  energy  ion  fluxes  are  strongly  anti-correlated 
with  9bx  and  essentially  little  correlation  for  energy  immediately  below  it.  It  is  even 
more  striking  that  this  critical  energy  is  identical  to  the  energy  at  the  spectral  break  of 
the  ion  spectra.  All  the  above  results  are  consistent  with  the  Fermi  acceleration  process 
since  it  is  responsible  for  the  energetic  ions  above  the  spectral  break  [e.g.,  Lee,  1982]. 
This  process  takes  place  at  the  quasi-parallel  bow  shock. 

Statistical  studies  show  that  probability  of  enhanced  energetic  ion  events  upstream 
and  downstream  from  the  bow  shock  increases  as  Obu  decreases,  and  vice  versa  [Mitchell 
and  Roelof,  1983;  Crooker  et  al,  1981].  In  principle  an  appropriate  Obu  should  be  used 
in  the  cross-correlation  analysis.  Since  the  observed  magnetosheath  energetic  ions  are 
extremely  energetic,  they  are  guided  more  or  less  by  the  magnetic  field  lines  not  the 
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flow  streamlines  in  the  subsonic  medium  [e.g.,  Bonifazi  and  Moreno,  1981;  Mitchell  and 
Roelof,  1983].  For  example,  a  41.1  keV/e  H'*'  or  He+^  ion  with  a  pitch  angle  of  45°,  the 
velocity  component  paxallel  to  the  magnetic  field  is  about  1900  km/s  which  is  more 
than  4  times  of  the  E  x  B  drift  speed  observed  in  the  magnetosheath.  Because  the 
associated  Osn  where  the  magnetosheath  magnetic  field  line  through  POLAR  intersets 
the  bow  shock  is  not  readily  known,  9bx  is  used  as  a  proxy  for  6 Bn-  Furthermore,  6bx 
and  the  angle  between  the  IMF  By  and  B^  components  would  completely  determine  the 
bow  shock  magnetic  geometry. 

Magnetosheath  energetic  ion  fluxes  are  enhanced  in  two  intervals,  ~9-10  and 
~11-12  UT.  Figure  10  shows  the  average  dayside  bow  shock  surface  projected  onto 
the  y-z  plane  for  the  first  interval.  This  bow  shock  surface  is  based  on  the  Fairfield 
model  [Fairfield,  1971].  For  the  IMF  observed  by  WIND  associated  with  this  interval, 
(—13.3,5.5,5.5)  nT,  6bx  is  30.3°  and  the  quasi-parallel  shock  as  shown  by  the  shaded 
region  covers  the  majority  portion  of  the  southern  hemisphere  and  part  of  the  northern 
hemisphere  of  the  dayside  bow  shock  surface.  During  the  second  interval,  IMF  Bx 
and  Bz  components  reverse  the  orientation  and  9bx  is  further  reduced  to  about  10°. 
The  quasi-parallel  bow  shock  switches  to  the  east  side  and  further  expands,  covering  a 
great  portion  of  the  southern  hemisphere  and  most  of  the  northern  hemisphere  of  the 
dayside  bow  shock  surface.  As  for  the  interval  of  weak  magnetosheath  energetic  ion 
fluxes,  ~10-11  UT,  9bx  is  about  80°  and  nearly  the  entire  dayside  bow  shock  surface  is 
quasi-perpendicular.  The  quasi-parallel  shock  in  this  case  occupys  only  a  small  portion 
of  the  dusk  flank  for  IMF  By  positive. 

The  majority  of  magnetosheath  energetic  ions  observed  by  POLAR  have  pitch 
angles  about  45°  (Figure  9).  Thus,  almost  all  the  energetic  ions  with  the  exception  of 
those  with  pitch  angles  nearly  90°,  would  tend  to  follow  the  magnetic  field  lines.  During 
the  first  interval  of  the  enhanced  magnetosheath  energetic  ion  flux,  ~9-10  UT,  both  IMF 
and  magnetosheath  magnetic  field  have  a  negative  Bx  and  positive  B^  component.  On 
the  basis  of  the  general  solar  wind  flow  and  magnetosheath  flow  [Spreiter  and  Stahara, 
1985],  the  most  likely  scenario  for  the  magnetic  field  topology  is  illustrated  in  Figure  11. 
POLAR  is  on  magnetosheath  magnetic  field  lines  that  are  connected  to  the  quasi-parallel 
bow  shock  (cf.  Figure  10).  Energetic  ions  produced  at  the  shock  follow  field  lines,  reach 
POLAR  and  are  detected  by  the  CAMMICE  instrument.  This  magnetic  topology  also 
appears  in  the  global  hybrid  simulations  results  (N.  Omidi,  private  communication).  For 
the  second  interval  of  intense  energetic  ion  flux,  magnetosheath  magnetic  field  Bx  and 
Bz  components  reverse  the  sign  but  the  magnetic  topology  remains  similar.  POLAR 
is  again  on  a  field  line  connected  to  the  quasi-parallel  bow  shock  and  observes  intense 
energetic  ions.  On  the  other  hand  between  the  two  intervals  of  enhanced  energetic  ion 
flux,  9bx  is  very  large  and  POLAR  is  on  field  lines  connected  to  the  quasi-perpendicular 
shock  and  detects  very  weak  energetic  ion  flux. 
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Figure  10.  Projection  of  the  dayside  bow  shock  surface  onto  the  y-z  plane  using  the 
Fairfield  model  [1971].  Contours  of  constant  Obt,  ‘‘•re  calculated  for  IMF  {B,,.,  By,  = 
(—13.3,5.5,5.5)  nT  with  9bx  =  30.3°.  The  shaded  region  is  the  quasi-parallel  shock 
located  at  the  most  part  of  the  southern  hemisphere  and  part  of  the  northern  hemisphere. 
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Figure  11.  Schematic  of  the  geospace  showing  the  bow  shock  (BS),  magnetosheath 
(MS),  magnetopause  (MP),  cusp  (CP),  and  the  POLAR  location.  POLAR  is  in  the 
magnetosheath  upstream  from  the  cusp  and  on  magnetic  field  lines  connected  to  the 
quasi-parallel  bow  shock,  observing  the  tailward-streaming  energetic  ions  that  are  accel¬ 
erated  at  the  quasi-parallel  shock.  The  magnetosheath  field  orientation  is  northward  and 
antisunward  in  favor  of  high-latitude  reconnection,  poleward  of  the  cusp  and  energetic 
magnetosheath  ions  can  then  directly  enter  the  cusp  along  the  interconnected  field  lines 
(e.g.,  field  line  c). 
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Because  the  observed  energetic  ions  have  a  very  high  speed  parallel  to  the  magnetic 
field  line,  it  takes  less  than  1  min  for  them  to  travel  from  the  quasi-parallel  bow 
shock  to  the  POLAR  location  in  the  magnetosheath.  The  solar  wind  propagation  time 
from  the  WIND  satellite  to  the  bow  shock  is  about  29  min  in  the  interval  selected 
for  the  correlation  analysis.  However,  the  associated  time  lag  for  the  best  correlation 
between  the  energetic  ion  flux  and  6bx  is  about  36-37  min  that  is  7-8  min  more  than 
the  above  solar  wind  propagation  time.  This  time  difference  can  be  attributed  to  the 
growth  time  for  the  energetic  ions  at  the  bow  shock.  As  bow  shock  changes  from 
the  cpiasi-perpendicular  to  the  quasi-parallel  geometry,  a  growth  time  of  ~10  min  is 
reciuired  for  enhanced  50-200  keV  ions  [e.g.,  Scholer  et  ai,  1980;  Mitchell  and  Roelof, 
1983].  Taking  into  account  the  uncertainty  of  a  few  minutes  (time  aliasing)  from  the 
CAMMICE  detector  sampling  time,  the  above  time  difference  is  consistent  with  the 
early  finding  of  the  average  growth  time. 

There  is  little  doubt  that  majority  of  magnetosheath  energetic  ions  came  from 
the  quasi-parallel  bow  shock  during  this  event.  Since  He"*"”  ions  are  exclusively  from 
the  Sun  not  the  ionosphere,  they  are  a  part  of  the  solar  wind  ions  accelerated  at 
the  quasi-parallel  shock  not  the  leaked  magnetospheric  ions  and  then  transported  to 
the  magnetosheath.  The  remarkable  similarity  between  the  magnetosheath  energetic 
H+  and  He+2  ions  in  the  spectral  shape,  correlation  coefficient  curves,  and  angular 
distribution  suggests  that  energetic  ions  are  likely  the  solar  wind  protons  accelerated 
at  the  cpiasi-parallel  shock  and  then  transported  to  the  magnetosheath  just  like  He'*'”. 
Magnetospheric  contribution  to  the  magnetosheath  energetic  ions  has  to  be  very  minor. 
The  simultaneously  enhanced  low-frequency  waves  associated  with  the  magnetosheath 
energetic  ions  (Plate  1)  are  likely  from  the  ciuasi-parallel  bow  shock  as  the  energetic 
ions  [e.g.,  Luhmann  et  ai,  1986;  Mdbius  et  ai,  1987].  These  magnetosheath  energetic 
ions  on  field  lines  close  to  the  bow  shock  eventually  cross  the  shock  and  escape  into 
the  interplanetary  medium  before  field  lines  are  connected  to  the  geomagnetic  field 
lines  following  magnetic  reconnection.  When  POLAR  was  upstream  from  the  cusp 
in  the  magnetosheath  and  magnetosheath  field  had  a  negative  B,;  and  positive  B., 
according  to  the  anti-parallel  merging  hypothesis  the  reconnection  site  was  probably  at 
the  high-latitude  magnetopause  poleward  of  the  cusp  [Crooker,  1979].  Magnetosheath 
field  lines  would  be  connected  to  the  cusp  field  line  and  energetic  ions  produced  at  the 
C[uasi-parallel  shock  would  then  enter  the  cusp  along  the  open  field  lines  (e.g.,  field  line 
c  in  Figure  11).  These  POLAR  observations  of  magnetosheath  energetic  ions  support 
the  hypothesis  of  bow  shock  source  of  CEP’s  [Chang  et  ai,  1998]. 
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7.  Summary  and  Conclusions 

We  have  presented  plasma  and  magnetic  field  observations  from  NASA 
GGS/POLAR  and  WIND  satellites  during  the  May  4,  1998  magnetic  storm  event. 
POLAR  was  at  the  dayside  magnetosheath  in  the  northern  hemisphere  and  observed 
energetic  (41.1-193.4  keV/e)  ion  fluxes  with  large  temporal  variations  (maximum  to 
minimum  flux  ratio  ~100).  Their  intensity  is  controlled  by  the  bow  shock  magnetic 
geometry.  These  ions  appear  to  be  from  the  bow  shock  accelerated  ions,  not  from  the 
magnetosphere  based  on  the  following  results. 

1.  POLAR  observed  distinct  ion  composition  (O'^''"^)  and  different  relative 
abundance  of  ions  (H"*",  He"^^,  and  0^"^^)  in  the  plasma  sheet  and  magnetosheath. 
Magnetosheath  energetic  ions  are  of  solar  wind  origin. 

2.  The  magnetosheath  ion  energy  spectra  are  continuous  with  a  spectral  break  at 
about  41.1  keV/e  for  all  three  ion  species,  resembling  those  at  the  quasi-parallel  bow 
shock.  Spectral  shape  for  energetic  and  He"^^  are  nearly  identical. 

3.  These  magnetosheath  energetic  ions  flow  antisunward  away  from  the  bow  shock 
along  the  magnetic  field  lines.  Angular  distributions  for  H"*"  and  He"*'^  are  very  similar. 

4.  Correlation  coefficient  for  the  magnetosheath  ion  flux  and  the  IMF  cone  angle 
(Obx)  depends  on  the  ion  energy  and  time  lag  for  For  energy  above  the  spectral 
break,  ion  flux  and  0£a:  offset  36-37  min  are  strongly  anti-correlated,  below  this  energy 
little  correlation. 

5.  Correlation  coefficient  curves  as  a  function  of  the  time  lag  peak  at  36-37  min. 
This  lag  is  consistent  with  the  solar  wind  propagation  time  estimated  by  four  different 
methods  plus  an  allowance  for  the  energization  time  of  the  energetic  ion  fluxes  at  the 
quasi-parallel  bow  shock. 

6.  Correlation  coefficient  curves  for  and  He"^^  with  energy  above  the  spectral 
break  are  nearly  identical. 

7.  When  POLAR  was  on  field  lines  connected  to  the  quasi-parallel  bow  shock, 
it  observed  intense  energetic  ion  fluxes.  When  it  was  on  field  lines  connected  to  the 
quasi-perpendicular  bow  shock,  httle  energetic  ions  were  observed. 

These  results  can  be  easily  understood  by  the  first-order  Fermi  acceleration  process 
at  the  quasi-parallel  bow  shock.  While  magnetosheath  magnetic  field  was  northward 
and  antisunward  in  favor  of  the  high-latitude  reconnection  poleward  of  the  cusp, 
energetic  magnetosheath  ions  upstream  from  the  cusp  were  streaming  tailward  in  the 
magnetic  field  direction.  As  the  field  lines  convected  and  were  interconnected  with  the 
geomagnetic  field  lines,  energetic  magnetosheath  ions  on  these  field  lines  would  travel 
along  the  interconnected  field  lines  into  the  cusp.  These  magnetosheath  energetic  ion 
observations  support  the  model  of  the  bow  shock  source  of  CEP’s  [Chang  et  ai,  1998]. 
Follow-up  studies  of  analyzing  multi-satellite  observations  at  the  bow  shock  and  in  the 
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cusp,  self-consistent  global  hybrid  simulations,  and  MHD  modeling  are  underway. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staffs  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analy¬ 
sis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid 
state  laser  design,  micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic 
frequency  standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation 
and  beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evalua¬ 
tion,  battery  electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  fail¬ 
ure  analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  space  environment 
effects  on  materials,  hardening  and  vulnerability  assessment;  contamination,  themal  and 
structural  control;  lubrication  and  surface  phenomena. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmos¬ 
pheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  sur¬ 
veillance,  imaging,  remote  sensing,  and  hyperspectral  imaging;  effects  of  solar  activity, 
magnetic  storms  and  nuclear  explosions  on  the  Earth's  atmosphere,  ionosphere  and  magne¬ 
tosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation,  design  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemdcal  reactions,  atmospheric  optics,  light  scattering,  state-specific  chemical 
reactions  and  radiative  signatures  of  missile  plumes. 

Center  for  Microtechnology:  Microelectromechanical  systems  (MEMS)  for  space 
applications;  assessment  of  microtechnology  space  applications;  laser  micromachining; 
laser-surface  physical  and  chemical  interactions;  micropropulsion;  micro-  and  nanosatel¬ 
lite  mission  analysis;  intelligent  microinstruments  for  monitoring  space  and  launch  sys¬ 
tem  environments. 

Office  of  Spectral  Applications:  Multispectral  and  hyperspectral  sensor  development; 
data  analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral 
imagery  to  defense,  civil  space,  commercial,  and  environmental  missions. 
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